The Antitriangular Factorization of Saddle Point Matrices

نویسندگان

  • Jennifer Pestana
  • Andrew J. Wathen
چکیده

Mastronardi and Van Dooren [this journal, 34 (2013) pp. 173–196] recently introduced the block antitriangular (“Batman”) decomposition for symmetric indefinite matrices. Here we show the simplification of this factorisation for saddle point matrices and demonstrate how it represents the common nullspace method. We show that rank-1 updates to the saddle point matrix can be easily incorporated into the factorisation and give bounds on the eigenvalues of matrices important in saddle point theory. We show the relation of this factorisation to constraint preconditioning and how it transforms but preserves the structure of block diagonal and block triangular preconditioners.

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Sparse block factorization of saddle point matrices

The factorization method presented in this paper takes advantage of the special structures and properties of saddle point matrices. A variant of Gaussian elimination equivalent to the Cholesky’s factorization is suggested and implemented for factorizing the saddle point matrices block-wise with small blocks of order 1 and 2. The Gaussian elimination applied to these small blocks on block level ...

متن کامل

Numerically stable LDL T - factorization of F - type saddle point matrices

We present a new algorithm that constructs a fill-reducing ordering for a special class of saddle point matrices: the F -matrices. This class contains the matrix occurring after discretization of the Stokes equation on a C-grid. The commonly used approach is to construct a fill-reducing ordering for the whole matrix and then change the ordering such that it becomes feasible. We propose to compu...

متن کامل

On the WZ Factorization of the Real and Integer Matrices

The textit{QIF}  (Quadrant Interlocking Factorization) method of Evans and Hatzopoulos solves linear equation systems using textit{WZ}  factorization. The  WZ factorization can be faster than the textit{LU} factorization  because,  it performs the simultaneous evaluation of two columns or two rows. Here, we present a  method for computing the real and integer textit{WZ} and  textit{ZW} factoriz...

متن کامل

THE USE OF SEMI INHERITED LU FACTORIZATION OF MATRICES IN INTERPOLATION OF DATA

The polynomial interpolation in one dimensional space R is an important method to approximate the functions. The Lagrange and Newton methods are two well known types of interpolations. In this work, we describe the semi inherited interpolation for approximating the values of a function. In this case, the interpolation matrix has the semi inherited LU factorization.

متن کامل

Null Space Algorithms for Solving Augmented Systems arising in the Mixed Finite Element Approximation of Saddle Point Problems

We use a Null Space algorithm approach to solve augmented systems produced by the mixed finite element approximation of saddle point problems. We compare the use of an orthogonal factorization technique for sparse matrices with the use of a sparse Gaussian elimination for the computation of the Null Space. Finally, we present the results of the numerical tests which we performed on selected pro...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

عنوان ژورنال:
  • SIAM J. Matrix Analysis Applications

دوره 35  شماره 

صفحات  -

تاریخ انتشار 2014